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Abstract —This paper deals with adaptive radar detection of 
a subspace signal competing with two sources of interference. 
The former is Gaussian with unknown covariance matrix and 
accounts for the joint presence of clutter plus thermal noise. The 
latter is structured as a subspace signal and models coherent 
pulsed jammers impinging on the radar antenna. The problem 
is solved via the Principle of Invariance which is based on 
the identification of a suitable group of transformations leaving 
the considered hypothesis testing problem invariant. A maximal 
invariant statistic, which completely characterizes the class of 
invariant decision rules and significantly compresses the original 
data domain, as well as its statistical characterization are deter¬ 
mined. Thus, the existence of the optimum invariant detector is 
addressed together with the design of practically implementable 
invariant decision rules. At the analysis stage, the performance of 
some receivers belonging to the new invariant class is established 
through the use of analytic expressions. 

Index Terms —Adaptive Radar Detection, Constant False 
Alarm Rate, Invariance, Maximal Invariants, Subspace Model, 
Coherent Interference. 

I. Introduction 

Daptive radar detection of targets embedded in Gaus¬ 
sian interference is an active research field, which has 
enjoyed the effort of world class scientists. As a result, 
in the last decades, a plethora of adaptive or nonadaptive 
detection strategies have been proposed to face with realistic 
applications. Generally speaking, the existing architectures can 
be classified into two macro-classes which differ in the level 
of a priori knowledge about the target response. The former 
assumes perfectly known target signature at the design stage 
(rank-1 receivers) Hl-lSl and references therein], whereas the 
latter takes into account possible uncertainties in the nominal 
steering vector Q-El. As a matter of fact, in situations 
of practical interest, several error sources (such as array 
calibration uncertainties, beampointing errors, multipath, etc.) 
ca can give rise to mismatched signals which can degrade 
the detection performances of the architectures belonging to 
the first class (7j, IfTTI . IflZl . 

The subspace detection represents an effective analytic tool 
incorporating the partial knowledge of the target response in 
the detector design and, hence, mitigating the performance 
degradation due to steering vector errors Ha-mu. The main 
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idea of the subspace approach is to constrain the target steering 
vector to belong to a suitable subspace of the observation 
space. By doing so, it is possible to capture the energy of 
the potentially distorted wavefront of a mainlobe target. As 
a result, architectures exploiting this trick are capable of 
declaring the presence of targets whose signature is signif¬ 
icantly different from the nominal one. It is worth pointing 
out that the subspace idea can be also exploited to model 
coherent interfering signals impinging on the radar antenna, 
whose directions of arrival have been estimated within some 
uncertainty. For instance, in ma, m the authors devise 
adaptive decision schemes to reveal extended targets when 
the interference comprises a random contribution (referring 
to clutter and thermal noise) and a structured unwanted com¬ 
ponent. 

In fI71 . 11221 . and fl4ll the subspace detection has been 
dealt with the invariance theory in hypothesis testing prob¬ 
lems l23l - |25l . By doing so, the authors focus on decision 
rules exhibiting some natural symmetries implying important 
practical properties such as the Constant False Alarm Rate 
(CFAR) behavior. Besides, the use of invariance leads to a 
data reduction because all invariant tests can be expressed in 
terms of a statistic, called maximal invariant, which organizes 
the original data into equivalence classes. Also the parameter 
space is usually compressed after reduction by invariance and 
the dependence on the original set of parameters becomes 
embodied into a maximal invariant in the parameter space (in¬ 
duced maximal invariant). Further examples of the application 
of such theory to radar detection problems can be also found 
in ma-Gi. 

It is now important to observe that some of the just 
mentioned papers do not assume the presence of structured 
interference IfTTI . f22|, l26l - lf33l . while others exploit the 
Principle of Invariance to solve the nonadaptive subspace 
detection problem when a coherent interference source illu¬ 
minates the radar antenna m. The possibility to address the 
adaptive subspace detection problem in the joint presence of 
random and subspace structured interference via invariance, 
to the authors best knowledge, has not yet been considered in 
open literature. In this respect, the scope of this paper is to fill 
this gap. To this end, this work is focused on adaptive detection 
of a subspace signal competing with two interference sources. 
The former is a completely random component, modeled as a 
Gaussian vector with unknown covariance matrix, and repre¬ 
sents the returns from clutter and thermal noise. The latter is a 
subspace structured signal (with unknown location parameters) 
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and accounts for the presence of (possible) multiple pulsed 
coherent jammers impinging on the radar antenna from some 
directions. The problem is analytically formulated as a binary 
hypothesis test and the Principle of Invariance is exploited at 
the design stage to concentrate the attention on radar detectors 
enjoying some desirable practical features. More specifically, 
a suitable group of transformation leaving the considered 
testing problem unaltered is established. Hence, the practical 
importance of the resulting group action is explained as a mean 
to impose the CFAR property with respect to the clutter plus 
noise covariance matrix and the jammer location parameters. 
Otherwise stated, all the receivers sharing invariance with 
respect to the devised group of transformations exhibit the 
mentioned CFAR features. Due to the primary importance of 
characterizing the family of invariant detectors, a maximal 
invariant statistic lf23l . lf24l . organizing the data into orbits 
and reducing the size of the observation space, is determined 
together with its statistical characterization. All the decision 
rules in the invariance family can be cast as functions of 
a maximal invariant. The problem of synthesizing the Most 
Powerful Invariant (MPI) detector l23l is also addressed and 
a discussion on the existence of the Uniformly MPI (UMPI) 
test is provided. Interestingly, the UMPI detector does not exist 
but for the case where the dimension of the useful signal plus 
jammer subspace equals the size of the observation space. In 
this specific situation, the energy detector turns out to be UMPI 
and coincides with the GLRT and the Locally MPI Detector 
(LMPID). When the UMPI test does not exist, the Conditional 
MPI Detector (CMPID), obtained as the optimum detector 
conditioned on the ancillary part of the maximal invariant 
is designed l23l . showing that it is statistically equivalent to 
the GLRT. At the analysis stage, we investigate the perfor¬ 
mances of some practically implementable invariant detectors 
providing analytical equations for the computation of the false 
alarm and the detection probabilities. The performance loss 
with respect to the clairvoyant MPI test is finally assessed 
showing that it is an acceptable value in some instances of 
practical interest. 

The paper is organized as follows. The next section contains 
the definitions used in the mathematical derivations, while 
Section III is devoted to problem formulation and derivation 
of the maximal invariant statistic. In Section IV, we provide 
the statistical distribution of the maximal invariant, whereas 
in Section V, invariant decision schemes are devised and 
assessed. Some concluding remarks and hints for future work 
are given in Section VI. Finally, Appendices lAllDl contain 
analytical derivations of the presented results. 


A. Notation 

In the sequel, vectors and matrices are denoted by boldface 
lower-case and upper-case letters, respectively. Symbols det(-) 
and Tr(-) denote the determinant and the trace of a square 
matrix, respectively. If A and B are scalars, then A x B is 
the usual product of scalars; on the other hand, if A and B 
are generic sets. Ax B denotes the Cartesian product of sets. 
The imaginary unit is j, i.e., y/—l = j. We denote by In 


the N x TV-dimensional identity matrix, while 0 is the null 
vector or matrix of proper dimensions. The Euclidean norm 
of a vector is denoted by || • ||. As to the numerical sets, IN 
is the set of natural numbers, 1R is the set of real numbers, 
R' VxM is the set of (TV x M)-dimensional real matrices (or 
vectors if M = 1), C is the set of complex numbers, and 
C ,v xM j s the set of (TV x M)-dimensional complex matrices 
(or vectors if M = 1). The subspace spanned by the columns 
of a matrix A £ C mxn is denoted by (A). The real and 
imaginary parts of a complex vector or scalar are denoted 
by Jf(-) and $$(•), respectively. Symbols (•)*, (-) T - and 
stand for complex conjugate, transpose, and conjugate trans¬ 
pose, respectively. The acronym iid means independent and 
identically distributed. If x £ C ,v x 1 is distributed according 
to the complex normal distribution with mean m £ C ,v x 1 and 
covariance matrix M £ (J NxN , we write x ~ CJ\f ^{m. M). 
The complex noncentral F-distribution with n, m complex de¬ 
grees of freedom and noncentrality parameter S is denoted by 
CJ- n ,m{5). Symbol CJ r „ im (0) is used to indicate the complex 
central F-distribution. The complex noncentral beta distribu¬ 
tion with n, m complex degrees of freedom and noncentrality 
parameter 5 is denoted by C/3 njm (<5), whereas we refer to 
the complex central beta distribution using symbol Cj3 n>rn ( 0). 
Symbol S ~ CWn(K, A) means that S £ <C NxN obeys the 
complex central Wishart distribution with parameters A' € IN 
and A £ C NxN positive definite matrix. Finally, if x £ € 
obeys the complex noncentral chi-square distribution with K 
complex degrees of freedom and noncentrality parameter S, 
we write x ~ Cx%{ 5); on the other hand, if x £ C is ruled by 
the complex central chi-square distribution with I\ complex 
degrees of freedom, we write x ~ Cxk( 0). 

II. Problem Formulation and Maximal Invariants 

In this section, we describe the detection problem at hand 
and introduce the Principle of Invariance that allows a com¬ 
pression of data dimensionality. Assume that a sensing system 
collects data from TV > 2 channels (spatial and/or temporal). 
The returns from the cell under test, after pre-processing, are 
properly sampled and organized to form a TV-dimensional 
vector, r say. We want to test whether or not r contains useful 
target echoes assuming the presence of an interfering signal. 

The target signature is modeled as a vector in a known 
subspace, (H) say, where H £ <C Nxr j r > 1, is a full- 
column-rank matrix. In other words, the useful echoes can 
be expressed as a linear combination of the columns of H, 
i.e., Hp with p £ C rxl . On the other hand, the interfer¬ 
ence component consists of two contributions. The former is 
representative of the clutter echoes and thermal noise, while 
the latter accounts for possible coherent sources impinging on 
the receive antenna from directions different to that where 
the radar system is steered. More precisely, the structured 
interferer signal is assumed to belong to a known subspace, 
(J) say, with J £ <C Nxt , t > 1, a full-column-rank matrix, 
and, hence, it is a linear combination, q £ C* x l say, of the 
columns of J with m = t + r < TV. The above model 
(with reference to Doppler processing) comes in handy to deal 
with scenarios where the presence of one or multiple coherent 
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pulsed jammers from standoff platforms attempt to protect a 
target located in the mainlobe of the radar antenna. We further 
assume that the matrix [J H] £ C ;V x m is full-column-rank, 
namely that the columns of J are linearly independent of those 
ofl H. Finally, as customary, we suppose that a set of K 
secondary data, r k £ <C Nxl , k = 1,..., K, K > N, namely 
data free of signal components and structured interference, is 
available. 

Summarizing, the decision problem at hand can be formu¬ 
lated in terms of the binary hypothesis test 


primary and secondary data can be transformed as follows 


z = Ur = U{[J H 


Q 

P 


n 0 


= U <Q [R q R p 


q 

p 


n 0 


— \Et E r \ (R q q + RpP ) + n 

9n 
02 

z k = Ur k =n k , k = l,...,K, 


— \Et E r 


+ n = E t Q i 


E r Q 2 


(5) 

(6) 


r = Jq + n Q , 

r k = n ok , k = 1,... ,K, 

(1) 

r = Hp + Jq + n 0 , 
r k =n ok , k=l,...,K, 

where 

• Hp £ C JVxl and Jq £ C ;Vx 1 are the target and interfer¬ 
ence signatures, respectively, with p and q deterministic 
and unknown vector^]; 

• no and not £ C Nxl , k = 1 are iid complex 

normal random vectors with zero mean and unknown 
positive definite covariance matrix M o £ <C NxN . 

Let us now recast i[7]i in canonical form lH3l . To this end, we 
exploit the QR-decomposition of the partitioned matrix [J H] 
given by 

[J H} = QR , (2) 



where Q £ <0 Nxm is a slice of unitary matrix (or a unitary 
matrix if m = N), i.e., Q^Q = I m (or Q^Q = QQ * = In 
if m = N ), and R £ C mxr " is a nonsingular upper triangular 
matrix. Moreover, observe that 


Q — [Q j Q 0 ]> n — 


Rj Rq 
0 Ri ’ 


(3) 


where Qj £ C Nxt and Rj £ € txt come from the QR- 
decomposition of J, namely, J = QjRj with Qj such 
that Q\Qj = It and Rj a nonsingular upper triangular 
matrix, Ro £ G txr , R\ £ C rxr ' is another nonsingular 
upper triangular matrix. Equalities are almost evident 
consequences of the Gram-Schmidt procedure. Now, let us 
define a unitary matrix U £ C NxN which rotates the columns 
of Q onto the first m elementary vectors of the standard basis 
of <C Nxl , i.e., 


UQ 


It 0 
0 I r 
0 0 


[E t E r \, 


(4) 


whcrcQ Et = [It 0 0] T and Et = [0 I r 0] T . Thus, under H\, 


is worth pointing out that the case t = 0 (i.e., absence of structured 
interference) has already been considered in open literature EJ. 

2 Observe that p also jointly accounts for target reflectivity, radar equation, 
and channel effects. 

3 It is worth noticing that when m = N, [Et E r \ = Ijy. 


where R q = [Rj 0] t , R p = [R q R-i] t , n = Un Q , n k = 
Un ok , k = 1,..., K, On = Rjq + R 0 p £ C txl , and 0 2 = 
Rip £ C rx 1 . On the other hand, when Ho is in force, the 
transformed data can be written as 


2 = Ur = UJq + n = E t Rjq + n = E t 6 10 + n , (7) 
z k = Ur k = n k , (8) 


where 0 io = Rjq £ C txl . Gathering the above results, we 
can rewrite problem ([Q in canonical form 


' H . f z = E t0 to + n, 

0 ’ \ z k = n k , k=l,...,K, 

< 

f z = E t 0 11 +E r 0 2 +n, 

1 ’ \ z k = n k , k = l,...,K, 

(9) 

where n, n k , k = 1...., K, are iid complex normal ran¬ 
dom vectors with zero mean and covariance matrix M = 
UM qU' 1 . Two important remarks are now in order. First, 
notice that (H) receives part of the interference energy and 
that (J) includes useful signal components since the columns 
H and those of J are not orthogonal but only linearly indepen¬ 
dent. This reciprocity becomes more evident after applying the 
transformations which lead to the canonical form. Second, in 
problem 0]) the relevant parameter to decide for the presence 
of a target is p. Otherwise stated, if the Hi hypothesis holds 
true, then ||p|| > 0, while ||p|| = 0 under the disturbance-only 
hypothesis (Hq). As a consequence, since f?i is nonsingular, 
problem © is equivalent to 


f H 0 : ||0 2 || = 0, 

1 Hi : || 0 2 || > 0 , 

which partitions the parameter space, © say, as 


( 10 ) 


0 = {0} U {0 2 g C rxl : ||e 2 || >0} . (11) 

00 ©1 

In the following, we can look for decision rules sharing 
some invariance with respect to those parameters (nuisance 
parameters, in this case M, On under Hi, and 0 io under IIq) 
which are irrelevant to the specific decision problem. To this 
end, we exploit the Principle of Invariance |}23l . whose main 
idea consists in finding transformations that properly cluster 
data without altering 

« the formal structure of the hypothesis testing problem 
given by H 0 : ||0 2 || = 0, JTi : ||0 2 || > 0; 

• the Gaussian model assumption under the two hypothe¬ 
ses; 
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• the subspace containing the useful signal components. 
Before introducing the transformation group that fulfills the 
above requirements, let us define Z s = [z± ... zk\. 



21 


' s 11 

Sl 2 

Sl 3 ' 

z = 

1 

CO to 

1 _ 

, 5 = Z S Z T S = 

1 

CO Co 

w to 

s 22 

S 32 

s 23 

S 33 


where Zl G C txl , z 2 £ C rxl , 23 G and S Kj , 

( i,j) G {1,2, 3} x {1, 2, 3}, is a submatrix whose dimensions 
can be obtained replacing 1, 2, and 3 with t, r, and N — m, 
respective! \Q. Moreover, Fisher-Neyman factorization theorem 
ensures that deciding from the raw data (z,Z s ) is tantamount 
to deciding from a sufficient statistic (z,S) (24) . 

Now, denote by Q£(N) the linear group of the N x N 
nonsingular matrices and introduce the following sets 


[ 

' G n 

g 12 

G 13 


\ G = 

0 

G 22 

g 23 

G G£(N) 

{ 

0 

0 

G 33 



G u G gc(t), G22 G gc(r), G3.3 G gc(N - m ) } (13) 


and 



/11 

0 


G C iVxl 


/n G € 



(14) 


It follows that the set Q x T together with the operation o 
defined by 


(G u fi) ° (G 2 , f 2 ) = (G 2 G 1 ,G 2 f 1 + / 2 ), (15) 

form a group, £ saj@, that leaves the hypothesis testing 
problem © invariant under the action l defined by 

l{z,S) = (Gz + f,GSG'), V(G, /) G £. (16) 

The proof of the above statement is given in Appendix [A] 
Moreover, it is important to point out that £ preserves the 
family of distributions and, at the same time, includes those 
transformations which are relevant from the practical point of 
view as they allow to claim the CFAR property (with respect 
to M and q) as a consequence of the invariance. 

Summarizing, we have identified a group of transformations 
£ which leaves unaltered the decision problem under consider¬ 
ation. As a consequence, it is reasonable to find decision rules 
that are invariant under £. Toward this goal, the Principle of 
Invariance is invoked because it allows to construct statistics 
that organize data into distinguishable equivalence classes. 
Such functions of data are called maximal invariant statistics 
and, given the group of transformations, every invariant test 
may be written as a function of a maximal invariant ESI. 


4 Hereafter, in the case m = N, the “3-components” are no longer present 
in the partitioning. 

5 Observe that 

• C is closed with respect to the operation o; 

. V(Gi, /r), (G 2 , / 2 ), and (G 3 , f 3 ) G £: [(Gi, JT) o (Ga, f 2 )\ ° 
(G? 3 , jT 3 ) = (Gi,/ 1 )o[(G, i / 2 )o{G 3 ,/ 3 )] (Associative Property); 

• there exists and it is a unique (Gy, f;) G C such that V(G, /) G £: 
(Gy, /y)o(G, f) = (G, /)°(Gy, fj) = (G, /) (Identity Element); 

• V(G, /) G C there exists (G_i, f _f) G £ such that (G_i, f _f) o 
(G, f) = (G, f) o (G_i, f_ 1 ) = (Gy, fj) (Inverse Element). 


The following proposition provides the expression of a 
maximal invariant for the problem at hand. To this end, it 
is necessary to distinguish between m < N and m = N\ the 
latter equality implies that the “3-components”, namely z 3 , 
S 23 , S 33 , S 31 , S 32 , G 13 , G 2 3 , and G 33 , are no longer 
present in the partitioned data matrices and vectors. 


Proposition 1. A maximal invariant statistic with respect to 
£ for problem © is given by 


ti(z,S) 


Zr > 3^2.3 Z 2.3 


"2.3 

-* 22 z 2 


Z 3 


m < N, 
m = N, 


(17) 


where z 2 .3 = z 2 -S 23 S 33 z 3 and S 2 .3 = S 22 -S 23 S 33 S 32 . 

Proof: See Appendix |B] ■ 

Some remarks are now in order. 

• Observe that when m < N the maximal invariant statistic 
is a 2 -dimensional vector where the second component 
represents an ancillary part, whose distribution does not 
depend on which hypothesis is in force. 

• It is important to compare (fTTb with the maximal invari¬ 
ant derived in El assuming the absence of subspace 
structured interference. Interestingly, they share a similar 
structure/dimensionality but for the presence of a projec¬ 
tion operation (onto the orthogonal complement of the 
rotated J) which amounts to cleaning the received data 
vectors from the structured interference. 

• Finally, note that (see Theorem 6.2.1 of |23| | ) every 
invariant test may be written as a function of (ITTI i (see 
Figure [T| for a schematic representation). 


III. Distribution of the Maximal Invariants 


In this section, we provide the statistical characterization of 
the maximal invariant for the case m < N and, then, we 
give a corollary referring to m = N. To this end, firstly 
we show that the maximal invariant can be written as a 
function of whitened random vectors and matrices and then we 
find a suitable stochastic representation by means of one-to- 
one transformations. Finally, the probability density functions 
(pdfs) in the invariant domain will be used in Section IV to 
construct the Likelihood Ratio Test (LRT) of Hq versus Hi, 
also referred to as the MPI Detector (MPID). 

Let us invoke the invariance principle and consider the 
following transformation (W, f 0 ) G £ 


where 


w = Vz + f 0 = 

S 0 = VSV t = 


V = 


V u 
0 
0 


Vi 2 

^ 22 

0 


w 1 




w 2 



(18) 

w 3 




Sou 

<S , 012 

S013 


St >21 

^022 

s 023 

, (19) 

£031 

S 032 

Sq33 



V 13 

—M 22 M 23 M 33 

m- 33 /2 


( 20 ) 
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In ©, V 22 = (M 22 - V ijt i = 1, 

j = 1,2,3, are matrices of proper dimensions, while M 22 £ 

<D rXr , M 23 S C rx(N-m)' and Ma3 g C (iV- m )x(tV-m) 

obtained partitioning M in the same way as S. 

Now, observe that under Hi, i = 0,1, 


w 2 

w 3 

S 022 S'o23 

§032 £033 


= V 2 


Z2 

Z3 


~ CN{iV 2 E r 0 2 ,I N -t) 
( 21 ) 



~CW N - t (K, 


S 23 
*S 3 3 . 

iN-t) 


V 


t 

2 


( 22 ) 


where 


V 2 


V 22 —M 22 M 23 M 33 

0 m - 3 1/2 


(23) 


Thus, we can recast ti(z, S) in terms of whitened quantities, 
more precisely, as 


m\ = {z 2 - S 23 S^z 3 )\S 22 - S 23 S^S 32 )~ 1 
x (z 2 — S 23 S 33 z 3 ) 

= ( wt 2 — *^ 023 *^ 033 ^ 3 )^ (*^022 — *^ 023 *^ 033 *^ 032 ) 1 

x (w 2 — Sq 23 S Q33 w 3 ), (24) 

m 2 = z\Sf 3 z 3 = ■mj5 ( 7 33 tn 3 . (25) 


Finally, apply a one-to-one transformation to (ITTb to obtain the 
following stochastic representation for the maximal invariant 


i 2 (z, S) = t 2 (w, S 0 ) 


1 


1 + 


m i 


1 + m 2 
1 


Pi 

P2 


1 + m 2 


(26) 

Summarizing, we have all the elements to provide the proba¬ 
bility density functions (pdfs) of the maximal invariant under 
both hypotheses. 


Proposition 2. The joint pdf of p\ and p 2 under the Hi 
hypothesis, i = 0,1, and assuming m < N has the following 
expression 


fpi,p^{x,y,Hi) = fp{x\K - (TV - f) + 1,r,0) 

x f p (y-K - (N - t) + r + 1, N - t - r, 0)e~ i SINR xy 

K — (N—t) + l 

x E 


fc =0 


K-(N-t) + 1 
k 


x (Xk-l )! [SINR (2?) 

where (x,y) £ {(0 1] x (0 1]}, 

SINR = d\(M 22 — M 23 Mf 3 M 32 )~ 1 6 2 (28) 


is the Signal-to-Noise-plus-Interference Ratio (SINR), and 
fp(x',n,m, 8 ) is the pdf of a l^j, / UJl/ 

• a random variable following a complex noncentral beta 
distribution with n,m complex degrees of freedom and 
noncentrality parameter 8 , if 8 > 0; 


• a random variable following a complex central beta 
distribution with n, m complex degrees of freedom, if 
8 = 0 . 

Proof: See Appendix [C] ■ 

Corollary 1. In the case m = N, the pdf of p 3 = 1/(1 + 
z\Sfl z 2 ) under the Hi hypothesis, i = 0,1, is given 


fp 3 i x > Hi) = fp{x;K -r + l,r,0)e 

K—r+1 


-i SINR x 


X 


E 

fc =0 


K — r + l\ (r — 1)! 


(r + k — 1)! 


[SINR (1 — x)] k , (29) 


where x £ (0 1] and SINR = 6 \M 2 2 0 2 . 


Proof: It is an application of the results contained in the 
proof of Proposition \2\ ■ 

As final remark, the induced maximal invariant can be easily 
obtained observing that the pdf of t 2 (z,S) depends on the 
unknown parameters only through the SINR. 


IV. Design and Analysis of Invariant Tests 

The goal of this section is twofold. First of all, we derive 
invariant decision schemes based on the LRT and its variants. 
In particular, we focus on the MPID, the LMPID (which is 
derived under the low SINR regime), and the Conditional 
MPID (CMPID), which are investigated for the case where 
the Uniformly MPID (UMPID) does not exist. Second, we 
consider some sub-optimum receivers providing the expres¬ 
sion of their decision statistics as functions of the maximal 
invariant. 


A. LRT-based Decision Schemes in the Invariant Domain 


In the sequel, we consider the cases m < N and m = N. 

1) Case m < N: As starting point of our analysis, we 
devise the MPID that, according to the Neyman-Pearson crite¬ 
rion, is given by the LRT after data compression. Specifically, 
the MPID has the following expression 


Lipid — 


fpi,P2 (•£> U i Hi ) 


fpitP2 (■**) v\ Hq) 


®=Pl.S(=P 2 


_ g—SINR P1P2 


(r-1)'. 

(r + k — 1)! 


fc=o 


K - (N - f) + 1 
k 


[SINR p 2 (l-pi)) h 


Hi 

J * 

H 0 


(30) 


where rj is the threshoIcQ set according to the desired Prob¬ 
ability of False Alarm ( Pf a ). It is important to observe 
that receiver is clairvoyant and, as a consequence, it 
is not of great practical interest in radar applications since 
it requires the knowledge of the induced maximal invariant. 
Nevertheless, it still represents a noteworthy benchmark, i.e., 
an upper bound to the detection performance achievable by 


s Observe that p 3 obeys the complex beta distribution with K r f I. r 
complex degrees of freedom. 

7 Hereafter, we use symbol rj to denote the generic detection threshold. 
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any invariant receiver. Before proceeding with the exposition, 
it is worth pointing out that decision scheme ( l30t is formally 
analogous to that found in G3, where the optimum detector 
is devised assuming a subspace structured signal embedded in 
random interference only (no subspace interference). This is 
perfectly consistent with the observation that when an additive 
structured interference is present, the MPID < l30b is obtained 
after projecting pre-processed data in canonical form onto the 
orthogonal complement of the interference subspace whose 
contribution is thus removed. 

Since in many applications the performances become criti¬ 
cal for low SINR values, the LMPID in the limit of zero SINR 
could be of interest. Following the lead of l34l and ll35l , the 
LMPID is given by 


t 


LMPID 


$fpi,P2(Pl,P2',Hl) 


(5SINR 

SINR=0 


fpi ,P2 (pi i Pz ■> Ho ) 


H i 

? * 

H 0 


(31) 


The following proposition provides the expression of f LMPID . 


Proposition 3. The LMPID is the following decision rule 


Karlin-Rubin Theorem, the test 

Ho 

Pi < r, (34) 

Hi 

is UMPI conditionally to p 2 . Moreover, we can consider the 
following equivalent form for test (l34l) 


t 


GLRT 


1 -Pi 
Pi 


Hi 

> 

< V, 

H 0 


(35) 


where, as we will show in the next subsection, f OLRT coincides 
with the GLRT in the original data space. 

2) Case m=N: Recall that when m = N the maximal 
invariant statistic obeys the complex beta distribution with 
K — r + l,r complex degrees of freedom and, hence, the 
LRT is given by 


fp 3 (*U HO 


K—r +1 

, __ J P3 v~> — -w __ -SINR p 3 

fp 3 K x 'i Hq) x=P3 

(r — 1)! . E 

X (^Wi)! [SINR(1 ^ 3)] < * 

Do 


K — r + 1 
k 


(36) 


tiMpm = —-—— ^ + 1 P2(1 - Pi) - PiPi < V- (32) 

H 0 

Proof: See Appendix O ■ 

As concluding part of this section, we provide the ex¬ 
pression of the CLRT, which is the LRT derived under the 
assumption that the ancillary part of the maximal invariant is 
assigned (23]. The main idea behind this approach relies on the 
fact that the ancillary part does not carry any information on 
the parameter and, hence, can be considered given. Moreover, 
it would be of interest to ascertain the existence of a CUMPID, 
since the UMPID does not exist for the problem at hand when 
m < N. 

The test statistic of the CMPID, t cMP id say, has the same 
expression as the MPID (see equation (l30t ) with the difference 
that in this case p 2 is a constant. Now, observe that i CMPID is 
a nonincreasing function of p\. As a matter of fact, let us 
evaluate those values of pi £ (0,1] such that the first derivative 
of f cmpid with respect to p\ is less than zero, namely 


St c 


Spi 


= —SINRp 2 e 

K-(N-t)+1 

x E 

k —0 


—SINR P 1 P 2 


K-(N-t) + l\ (j— 1)! 
k 


E 

fc =i 


(r + k — 1)! 

SINR P 1 P 2 

K-(N-t) + l\ (r — 1)! 


x [SINR p 2 {l-pi)} k - SINRp 2 e" SINR PlP2 

K-(N-t) +1 


(r + k — 1)! 


x [SINR p 2 (l — Pi)]- 1 < 0. 


(33) 


The above inequality holds true Vpi £ (0,1] and, hence, f C MPiD 
is a nonincreasing function of pi. As a consequence of the 


Now, following the same steps used to prove that the GLRT 
is conditionally UMPI, it is easy to show that the test 

Ho 

P3 < V (37) 

Hi 

is UMPI and statistically equivalent to the Energy Detector 
(ED) given by 

1_ Hi 

tm= -- = z\ S 22 z 2 < V- (38) 

P3 A 

Do 

Finally, the following corollary provides the expression of the 
LMPID when m = N. 

Corollary 2. The expression of the LMPID in the case m = N 
is 

K — r + 1 

-(1 P 3 ) P3 < V- (39) 

H 0 

Proof: The above result can be easily obtained following 
the line of reasoning of the proof of Proposition [3 ■ 

It is clear that the LMPID for m = N is statistically 
equivalent to the ED. 

B. Receivers Obtained as Functions of the Maximal Invariant 

Several invariant decision rules can be obtained exploiting 
different rationales based on either solid theoretical paradigms 
or heuristic criteria. The basic observation is the property that 
any invariant test can be written as a function of the maximal 
invariant. Thus, it is possible to device invariant architectures 
choosing proper functions of the maximal invariant. Specifi¬ 
cally, we focus on the GLRT, which under some mild technical 
conditions always leads to an invariant detector 1361 , and the 
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so-called two-step GLRT (2S-GLRT) J2] which consists in 
devising a non-adaptive GLRT for known covariance matrix 
(step one) and, then, in replacing the unknown covariance 
matrix with a proper estimate (step two) to come with a fully 
adaptive decision scheme. Again, in the sequel, we distinguish 
between the cases to < TV and to = TV. 

1) Case m < TV; Let us begin by providing the test statistic 
of the GLRT, which has the following expression E9, ED 


l + ziS-^Pj-y^S-^z 

1 + zt S~ l/2 Pg_ l/2Em S~ l/2 Z 


(40) 


where P 4 = A(A 1 A)^ 1 A' 1 is the projection matrix onto the 
subspace spanned by the columns of A £ (C ,VxM , P 4 = 
In — Pa, and E m = [ E t E r ], It is tedious but not difficult 
to show that 

z'S- 1/2 Pj. 1/2E S~ 1/2 z = 1 ~ PlP2 , (41) 

* P 1 P 2 

z^S~ 1/2 Pg- 1/2E S~ 1/2 z = (42) 

p 2 

Sumarizing the statistic of the GLRT can be recast as 

foLRT = ) (43) 

Pi 

which is clearly statistically equivalent to the left-hand side of 

(E3. 

On the other hand, the two-step GLRT (2S-GLRT) is given 
by 

W rt = z^S~ 1/2 Pj. 1/2Et S~ 1/2 z 

— 2 , S” 1/2 PL / 2 £ S~ 1/2 z = (44) 

P 1 P 2 

It follows that both receivers are invariant with respect to C 
and, hence, they ensure the CFAR property with respect to the 
unknown covariance matrix of the disturbance and q. 

2) Case m = TV; In this case, the GLRT and the 2S-GLRT 
share the same test statistic which has the following expression 


fr.i.RT — t 


GLRT — L 2S-GLRT 


— z t C-t/2 p-L Cf-1/2 


Z. 


It is not difficult to show that 

>-L 

S~ 1 / 2 E, 


te-!/2p-l- q-1/2 _ t c-t. _ 

Z a * <Z- 1 / 2 B.O Z — Z 2 J 22 Z 2 — 1 

P3 


(45) 


(46) 


which is the statistic of the ED. It follows that, under the 
assumption m = N, the GLRT, the 2S-GLRT, and the LMPID 
coincide with the ED, which, in turn, is UMPI. 


1) Case m < N: Let us focus on the GLRT and observe 
that under the H 0 hypothesis, f GLRX ~ CJ- r x-(N-t)+ 1(0), 
while under Hi, t 0LRT ~ CP rK _^ N _ t ) + i(S) with S 2 = 
P 2 SINR. It follows that the Pf a can be expressed as 


P?T(V) = 1 - F(q; r ,K-(N-t) + 1, 0) 

r —1 


1 


(1 + v ) r + K - N + t 


E 

;=o 


r + K — N + t 
l 


(47) 


where F(x\n,m,8), x > 0, is the Cumulative Distribution 
Function (CDF) of a random variable ruled by the complex 
central (noncentral) F-distribution with n, m complex degrees 
of freedom if S = 0 (<5 > 0). Following the same line of 
reasoning as for the Pf a , it is not difficult to show that /-)/ is 
given by 


Pr T (v, SINR) = 1 - ^ F( m r ,K-(N-t) + 1, S(u)) 
Jo 

x fp(w, K — (TV — t) + r + 1, TV — t — r, 0 )du, (48) 


where <5 2 (it) = uSINR. 

As to the 2S-GLRT, observe that it can be recast as Ls-glrt = 
Tglrt/ P‘ 2 , and, hence, Pf a and Pd can be easily derived relying 
on the above results, more precisely 


Pf S a LRT (v) = 1 - T F(vu-, r, K - (TV - t) + 1, 0 ) 

Jo 

x fp{u\ K — (TV — f) + r + 1, TV — t — r, 0 )du, 

(49) 


Pr LRT (P, SINR) = 1- f 1 F^ryu-, r, K — (N — t) + l, 5(u)) 
Jo 

x fp(u\K — (TV — t) + r + 1, N — t — r, 0 )du. (50) 


In order to provide the expressions of Pf a and Pd for the 
C. Pfa and Pd of the GLRT, the 2S-GLRT, the LMPID, and LMPID, let a = \K - (N — T) + 1 ]/r and observe that 
ED 


This subsection is aimed at the derivation of closed-form 
expressions for the Pf a and the Probability of detection (Pd) 
of the GLRT, the 2S-GLRT, and the LMPID. To this end, we 
deal with equivalent decision statistics and distinguish between 

to < TV and to = TV. 


ap 2 {l - Pi) - P 1 P 2 = P 1 P 2 



P1P2 [aWT - 1] = 


P2 


^GLRT + 1 


( Uf GLRT 1) ■ 



(51) 











This implies that Pf a can be written a 


HI 


distribution of f OLRT , namely 


PfT D (v) 

Vn 

= P J (ctfoLRT 1 ) > it, Ho 

Tglr'i T 1 

= p [fGLRx(ap 2 - v) > P 2 + rn H 0 ] 

= / P [toLRiiau — rj) > u + rj\p 2 = u; H 0 \ 

Jo 

x fp(w, K — (TV — t) + r + 1, N — t — r, 0 )dv 

' f 1 r_ u + r] 

/ P ioLRT > - \P 2 = u; H 0 

Jri/a L au — rj 


= < 


X 


+ / p 

Jo 


— XL + TJ 

Tglrt < P2 = U", Ho 

au — p 


(52) 


rj/a L uu- 1/ 

fp{u\ K — (N — t) + r + 1, N — t — r, 0)du 

rri/a r . 


- r + 1, TV — f — r, 0)du, 


J o |_ CLLL f j 

x/^( it; K — (TV — t) + r + 1, TV — f — r, 0)dzi, 
if 0 < 77 < a, 

0 , if 77 > a, 

f ^ ^ _ j_ ^ 

/ P Tglrt > |P2 = TTo 

jo L au - V 

xfp^K - (N — t) +', 

if — 1 < r] < 0 , 

1 , if 77 < — 1 , 

f \l - F (^--r,K - (N-t)+ l,o) 
J V /a l \Ua-T) J 

xfp(u; K — (TV — f) + r + 1, TV — t — r, 0 )du 

+ [ V/a F(^L;r, K-(N-t)+ l,o) 

Jo \ua-r) J 

xfp(u ; TV — (TV — t) + r + 1, TV — t — r, 0)d«, 
if 0 < 77 < a, 

0 , if 77 > a 

1- I' F(^L;r,K-(N-t)+ 1,0) 

Jo \ua — rj J 

xfp(u ; K — (TV — t) + r + 1, TV — t — r, 0)du, 
if — 1 < 77 < 0 , 

1 , if 77 < — 1 . 


The expression of the /-)/ is formally analogous to that of the 
Pf a but for the presence of the noncentrality parameter in the 


8 Notice that — 1 < tunpro < a - 


P^ mpid (t 7 , SINR) 
rl 

lr)/ 


= 


It + 77 


; r, K — (TV — t) + 1,8 {u) 


ua — 77 

xfp(u; K — (TV — f) + r + 1, TV — t — r, 0)du 


rv/a 


F 


u + 77 


; r, K — (TV — t) + 1, 8 (u) 

ua — 77 

xfp(u] K — (TV — f) + r + 1, TV — t — r, 0)du 
if 0 < 77 < a, 

0 , if 77 > a, 

r 1 / , 1 

U -;r,K — (N — t)+ l,S(u) 


1 - 


F 


v zta — 77 

xfp{u\ K — (TV — f) + r + 1, TV — t — r, 0)d7t 
if — 1 < ?7 < 0 , 

1 , if 77 < —1. 


(54) 


2) Case m = TV; Under the assumption m = TV, the GLRT, 
the 2S-GLRT, and the LMPID are statistically equivalent to the 
ED. For this reason, the object of this section is the statistical 
characterization of the ED. Specifically, it is easy to show that 
tuD is ruled by 

• the complex central F-distribution with r, K — r + 1 
complex degrees of freedom under Ho', 

• the complex noncentral F-distribution with r,K — r + 1, 

degrees of freedom and noncentrality parameter 5 with 
S 2 = under H\. 

It follows that 



and 

PTiV, SINR) = 1 - F( 77 ; r,K-r + l,S). (56) 


D. Illustrative Examples 

Figures [2]|5] display /)/ versus the SINR for the considered 
decision schemes. The curves have been obtained by means 
the above closed-form formulas which have been computed 
by standard numerical routines but for the MPID, whose 
performance has been obtained via standard Monte Carlo 
counting techniques. More precisely, the thresholds necessary 
to ensure a preassigned value of Pf a have been evaluated 
exploiting 100/P/ a independent trials, while the Pd values are 
estimated over 5000 independent trials. As to the disturbance, 
it is modeled as an exponentially-correlated Gaussian vector 
with covariance matrix M = aalN+<J 2 M c , where > 0 is 
the thermal noise power, a 2 > 0 is the clutter power, and the 
(i, j)-th element of M c is given by 0.95l* _J L The clutter-to- 
noise ratio cr^/cr^ and the interferer-to-noise ratio || 0 i|| 2 /cr^ 
are both set to 30 dB with a 2 = 1. Finally, all numerical 
examples assume P/ a = 10~ 4 . 

Figures [2] and [3] show the performances of the considered 
detectors assuming TV = 8 , K = 12, and different values 
of (r, £). Inspection of the figures highlights that the GFRT, 
the 2S-GFRT, and the MPID share practically the same 
performance, while the FMPID performs poorer. In particular. 
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the latter exhibits a loss of about 3 dB at Pd = 0.9 in 
Figure U and does not achieve Pd = 0.9 for the SNR values 
considered in Figure [3] Moreover, the GLRT and the MPID are 
slightly superior to the 2S-GLRT. Finally, it is worth observing 
that for high values of r, the Pd curves of the considered 
detectors move toward high values of SINR, namely there 
is a performance degradation with respect to the case where 
r is low; this is strictly tied to the increase of the number 
of unknown parameters to be estimated. In Figures Q] and [3 
Pd versus SINR is plotted assuming N = 8 and K = 16. 
Inspection of the figures highlights that the above hierarchy 
keeps unaltered, while the increased number of secondary data 
used for estimation purposes decreases the loss of the 2S- 
GLRT with respect to the GLRT. Finally, notice that contrary 
to the previous case the Pd of the LMPID can be higher than 
0.9 for SINR values greater than 16 dB when r = 2 and 18 
dB when r = 4. 


V. Conclusions 

In this paper we have considered adaptive radar detection 
of a subspace signal embedded in two sources of interfer¬ 
ence: a Gaussian component with unknown covariance matrix 
modeling the joint effect of clutter plus receiver noise and a 
subspace structured part accounting for the effect of coherent 
jammers impinging on the radar antenna. We have formulated 
the problem as a binary hypothesis test and have exploited the 
theory of invariance to characterize adaptive detectors enjoying 
some relevant symmetries. In this context the main achieved 
technical results are: 

• The identification of a suitable group of transformations 
leaving the considered hypothesis test invariant and the 
interpretation of the group action as a mean to force at 
the design stage the CFAR property with respect to clutter 
plus noise covariance matrix and jammer parameters. 

• The design and the statistical characterization of a max¬ 
imal invariant statistic which organizes the radar data in 
equivalence classes and significantly realizes a compres¬ 
sion of the original observation domain. 

• The synthesis of the optimum invariant detector and the 
discussion on the existence of the UMPI test. In this 
respect, we have shown that the UMPI detector in general 
does not exists but for the case where the useful signal 
plus jammer subspace size completely fills the dimension 
of the observation domain. 

« The synthesis of the LMPID, the CUMPID, and the proof 
that the latter test coincides with the GLRT. 

• The development of analytic expressions for the per¬ 
formance assessment of some practically implementable 
invariant detectors. 

Possible future research tracks might concern the possibility 
to deal with a partially homogeneous Gaussian interference 
through the additional invariance under a common scaling of 
the secondary data, as well as the eventuality to force some 
structure in the covariance matrix of the Gaussian interference 
(as for instance persymmetry (37), (38l ). Last but not least, 
it would be interesting the extension of the entire framework 
to a non-homogeneous scenario where the training data may 
exhibit different power levels. 


Appendix A 

Invariance of Problem ® with respect to the 
Group C 


Let (G, /) € C and observe that, under Hi, x = Gz + f 
obeys the complex normal distribution with covariance matrix 
GMG t and mean 


G(E t e u 


/11 




' Guflu ' 


G1202 

E r e 2 ) + / = 

0 

+ 

G 22 G 2 


0 


0 


— EfO i + E r 0 2 , (57) 


where 6 \ = GnOn + G 12^2 + fn and 62 = G 22 $ 2 - 
Moreover, x is independent of x 1 = Gzk, k = 1..... K, 
which are iid complex normal vectors with zero mean and 
covariance matrix GMG Now, observe that |6L|| > 0 if 
and only if ||^ 21 | > 0 (since G 22 is nonsingular). 

On the other hand, when Hq is in force, x is distributed as 
under H 1 but for the mean, which is given by 


GE t e w + f 


GnOio + fn 

0 

0 


E t G 10 , (58) 


where #io = GnO 10 + / n , whereas x &, k = 1,..., K, are 
still iid complex normal random vectors with zero mean and 
covariance matrix GMG^. 

Thus, we have shown that the original partition of the 
parameter space, data distribution, and the structure of the sub¬ 
space containing the useful signal components are preserved 
after the trasformation (G, /). 


Appendix B 

Derivation of the Maximal Invariant Statistics 

This appendix is devoted to the proof of Proposition [7] In 
particular, we focus on the more challenging case 2 < m < N 
since m = N can be obtained according to the same proof 
method. 

Before proceeding, it is worth recalling that a statistic 
t(z,S) is said to be maximal invariant statistic with respect 
to a group of transformations C if and only if (23l . (251 

t(z,S)=t[l(z,S)\, VZ G £, (59) 

and 


t(z,S)=t(z,S)=>3l&£ : (z,S)=l(z,S). 

In order to prove (f59l) . we consider the following vector 

z 23 ^ 2^ z 23 

. Z 3 <S 33 z 3 . 


(60) 


(61) 


where 


223 = 


22 

23 


, s 2 = 


S 22 S 23 

S 32 S 33 


and 


z 23 S 2 1 2'23 = Z \S^Z 3 + Z2.3S2lz2.3- 


(62) 

(63) 
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Moreover, let us partition the matrix G as follows 
G = 


G\ G 2 

0 g 3 


G, = 


G 22 G 23 

0 G 33 


G 1 Z 1 + G 2 Z 23 + f 11 
G 3 Z 23 


Now, let (z, S) = l(z , S) and observe thaf| 
z = Gz + f = 

and 

S = GSG ’ = g 3S ,gU 

where Ai, i = 1,2, 3, are matrices of proper size and 


g 3 s 2 g\ = 


B 1 


B, 


B 3 G 33 S 3 3 G 33 


t 


z\ 3 S 2 z 2 3 

^3^33 z 3 

Z 23*^2 Z 23 
. z 3 . 


z\ 3 Gl{G\)- 1 S 2 1 G 3 1 G 3 Z23 

z\G33(G1 3 )- 1 S 3 ^G^G33Z3 


and, hence, that 


^2.3^2.3^2.3 — Z 2.3'S'2.3 Z 2.3- 


Z 2.3^2.3 Z 2.3 


z 2.3^2.3 z 2.3 

z 3'^33 2 3 


z 3^33 z 3 


I2/2.3II 2 = fej 2 , 

II 2 / 3 II 2 = lly 3 ll 2 , 


Moreover, partition S 2 1 and S 2 1 as follows 


(64) 


SP = 


2 a,,u 2 

I r 0 

33 S 32 I N—m 


* 2 .: 

0 


S23 0 

S 33 1 


where Gi = G n G C txt , G 2 = [G 12 G 13 } G C* x and 


G C (jv-t)x(tv-t)_ (65) 


P f 

Ir -S 2 3 S 3 i 

0 IN—m 


S 2 1 = 


P 

I r 0 

33 J 32 iN—m 


Son S, 


w 2 


5 1 

5 2.3 

0 S' 


(75) 


0 

-1 

33 


( 66 ) 


(67) 


( 68 ) 


>t 


Ir -S 23 S , 1 


w 2 


’23^33 
IN—m 


(76) 


and note that 


with, in turn, B t , i = 1,2,3 matrices of proper dimensions. 
It follows that 


'23 


'23 



C — 1/2 r 
‘-*2 3 Z 2.3 


2/2.3 


Q^ 1 / 2 - 
*33 Z 3 


2/3 


^/ 2 y 

*2.3 z 2.3 


2/2.3 


o^l/2 - 
*33 z 3 


Va 


Gathering the above results yields 
WPZ 23 = 


(69) 


(70) 


U 2.3 0 

0 u 3 


WPz 


23 


u l 

Z23 = (wpy'ihwP^ 


(77) 

(78) 

(79) 

(80) 


The proof of the second property (fhOl i adopts the following 
rationale: first, we find the submatrix G 3 and then we construct 
/ n and the remaining blocks of G. 

Assume that there exist (z, S ) and (z, S) such that 


where D 3 = (WP) _ 1 £/ i\V P is an upper block-triangular 
matrix with the same structure as G 3 in ( 165 b . Besides, from 
equations (1751 ) and (1761) it follows that 


WPS 2 P ] W = I N —m = UiU[ 

=> U\WPS 2 (U\WPy = I N —m = WPS 2 {Wpy 


t 


(81) 


t( z , s) = = 2 ?‘=-?3 = t (z,s). 

(71) 

The above equalities can be recast using the Euclidean norm 
of a vector, more precisely 


(72) 

(73) 


(82) 

=► S 2 = ( WP)- 1 U 1 WPS 2 {Wpyu\[{WPy ]- 1 (83) 
= D 3 S 2 DI (84) 

So far, we have constructed the block G 3 = D 3 of G, it 
still remains to find the other blocks of G. To this end, let 
S3 = [S 12 S' 13 ], S 1 = S 11 and write d67l > as 


G~ 1 S = SG i 


(85) 


where y 2 3 = S 2 3 7 z 2 . 3 , V 2 .3 = S 2.3 z 2 . 3 , 2/3 = S 33 7 z 3 , 
and y 3 = S 33 z 3 . As a consequence, there exist U 3 £ 
(^(jv-m)x(jv-m) an( j JJ 2 3 g C rxr ' un j tar y matrices such that 

y 2 .3 = U 2 .3V2.3 and 2/3 = G 3 y 3 . ( 74 ) 

9 We partition z and S according to the same rule used for z and S. For this 
reason in the sequel we omit the definition of the submatrices and subvectors 
of z and S. 


G^ 1 ~G 3 1 G 2 G 3 1 

0 G7 1 


' S 1 S 3 


' s t S3 ■ 


' g\ 0 

s\ s 2 


_ S 3 s 2 


. G't G 3 _ 


( 86 ) 

G^Si - Gi 1 G 2 G 3 1 Sl Gf 1 ^ - Gi L G 2 G 3 1 S 2 

g 3 1 s 2 


G^Sl 


S'lGl + S 3 G\ S 3 Gl 

sl g\ + s 2 g\ s 2 g\ 


( 87 ) 
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The last equation is equivalent to the following system of 
equations 

{ G7 1 S 1 -G7 1 G 2 G7 1 S\ = S 1 G] + S 3 gI, 

| G/ 1 ^ - G^ 1 G 2 G^ 1 S 2 = S 3 g\, (88) 

[ = S'gGj + S 2 G\. 

Replacing the following equal inf'] 

G^s 2 = S 2 G\ (89) 


into the second equation of the above system yields 


G 3 l S 3 - G^G 2 S 2 G\ = S 3 g\ (90) 

=* S 3 - G 2 S 2 gJ = GrSaGj (91) 

^(G^ 1 ^) t = G 1 5 3 + G 2 5 2 (92) 

=► G^Sl = 5^Gj + S 2 G\. (93) 


Thus, the second matrix equation of (188b is redundant and can 
be neglected. From the third equation it stems that 


g\ = s ; 


1 (gj 1 ^ - S^Gj 


which, replaced in the first equation of 


leads to 


(94) 


G^Si-Gr 1 S 3 (Gl)~ l - G 3 S 3 


S^G^Sl 


= SiG\ + S 3 S 2 \G 3 1 Sl - S' 3 G\ 

=> S t - SsS^Sl + G^sS^G^Sl 




(95) 


= GiS'iGj + G 1 S 3 S 2 1 GJ 1 SI - G 3 S 3 S 2 ^Gj (96) 
=► Gi(Si - SsS^S^Gl = Si - SsS^Sl (97) 


The solution of the last equation has the following expression 

Gi = (Si - S 3 S 2 1 SJ) 1/2 (Si - SsS^sl)- 1 / 2 (98) 


and it is such that det(Gi) ^ 0. Replacing Gi in (|94] > with 
( |98| > yields the expression of G 2 . 

Finally, / n can be evaluated as follows 


f — z — Gz = 


Zi 


Gi 

g 2 


Z 1 

Z 23 


0 

g 3 


z 23 _ 


z\ — G\Z\ — G 2 z 23 

0 


(99) 

( 100 ) 


where the last equality comes from ( l80l >. 


Appendix C 

Statistical Characterization of the Maximal 
Invariant 

The joint pdf of j> \ and p 2 can be computed exploiting the 
following equality 

fpi, P 2 (x,y) = f P 1 \ P 2 ( x \p 2 = y)f P 2 (y), ( 101 ) 


Let us focus on f P2 ( y ) and observe that from (OH it stems 
that w 3 ~ CAf N -m( 0 , I N - m ), while by Theorem A. 11 of 0 
S 033 ~ CW N -m(K,I N - m ). Now, recast w\Sq 33 w 3 as 


P 2 


wlw 3 

wj W 3 

W 3 Sq 33 W 3 


a 

b 


( 102 ) 


where a ~ Cx^_ m (0) and, by Theorem A.13 of 0 , b ~ 
CXK-(N-m)+ 1(0) is statistically independent of a. It follows 
that p 2 is ruled by the complex central F-distribution with 
TV — to, K — (TV — m) + 1 complex degrees of freedom and, 
hence, p 2 = 1/(1 + p 2 ) ~ C/3 K _( N _ m ) + i,;v_ m (0). Note that 
p 2 is an ancillary statistic and its distribution is the one and 
the same under both hypotheses. 

The distribution of p 3 can be obtained assuming that the 
“3-components”, namely the entries of the considered random 
vectors sharing the subscript equal to 3, are no longer random 
variables but are assigned. In such a case, p 2 is deterministic 
and under Hi 

. d = (w 2 - s023 Sq 33 W 3 )/^J 1 + wISm 3 w 3 ~ 

CU r ((M 22 - M 23 M33 1 M 32 )- 1 / 2 0 2l / r ); 

. X = (S'o 2 2 - S 023 So 3 1 3 S 032 ) ~ CW r {K - (TV - t) + 
r, I r ) independent of d (see Theorem A. 11 of 0 ). 
Again, exploiting Theorem A.13 of 0 leads to 

Pi =-Tjr-— ~ C.7F r ,if-(jv-t)+i (^)i 

d'd 

d)X~ x d 

where 

^ 2 __ 0 2 (M 22 — M 23 M 3 ^M 32 ) 1 0 2 

l + wls 033^3 


(103) 


(104) 


is the noncentrality parameter. As a consequence, given p 2 . 
Pi = 1/(1 +pi) ~ C(3 K _ {N _ t)+hr (6) with b 2 = SINRp 2 . It 
is clear that under H 0 and given p 2 , p 3 ~ C/3x-(A-t)+i,r(0), 
which does not depend on p 2 . 

Gathering the above results, the joint pdf of p 3 and p 2 can 
be written as 


f P i, P 2 ( x ,y,Hi) = fp{x\K- (TV - t) + 1,r, 0) 
x fp{y, K ~ (TV — f) + r + 1, TV — f — r, 0) 

K-(N-t)+l 

x e _i SINR xy 

fe=o 

(r-1)! 


K - (TV - t) + 1 
k 


(r + k — 1)! 


[SINR y{ 1 - x)}\ (105) 


under 11 \, and 


fpi,p 2 { x ,y,Ho) = fp(x;K- (TV — f) + 1, r, 0) 

x fp{y ; K — (TV — t) + r + 1, TV — t — r, 0), 

(106) 


where f Pl \ P 2 (x\p 2 = y) is the conditional pdf of p 3 given p 2 
and f P 2 {y) is the marginal pdf of p 2 . 

10 Observe that it comes from EJ. 


under Hq, where 

fp{x; n, ro, 0) = + x ^ 1 {l-x) m ~ l . (107) 

(n — l)!(m — 1 )! 
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Appendix D 

Derivation of the LMPID 
It is not difficult to show that 


K-(N-t) +1 
a -SINR pip2 y' 

fc =0 

[SINR P2 (l-Pi)] k 


) (5SINR 

(r~ 1 )! 

(r + k — 1)! 

( K-(N-t) + l 

— piP2e~ slNR PlP2 

l fe=0 


K-(N-t) + 1 
k 


SINR=0 

K - (N - f) + 1 
k 


x ( r . + fc i ) i)| t SINR ^(l-Pi)] fc +e 


—SINR P 1 P 2 


iC-(AT-£) + l 

x E 

fc=i 


if - (iV - t) + 1\ (r — 1)! 

k J (r + k — 1)! 


x fcSINR fc_1 [p 2 {l-pi)} k 
K - (N - t) + 1 


) SINR=0 
P2(l -Pi) -PlP2- 


(108) 
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Fig. 1. Block diagram of the transformations that lead to a generic invariant 
test. 



SINR (dB) 

Fig. 2. Pd versus SINR for the GLRT, the 2S-GLRT, the LMPID, and the 
MPID assuming N = 8, K = 12, r = 2, t = 4, and Pf a = 10 -4 . 


Fig. 3. Pd versus SINR for the GLRT, the 2S-GLRT, the LMPID, and the 
MPID assuming N = 8, K = 12, r = 4, t = 2, and Pf a = 10 -4 . 



Fig. 4. Pd versus SINR for the GLRT, the 2S-GLRT, the LMPID, and the 
MPID assuming N = 8, K = 16, r = 2, t = 4, and Pf a = 10 -4 . 



SINR (dB) 


Fig. 5. Pd versus SINR for the GLRT, the 2S-GLRT, the LMPID, and the 
MPID assuming N = 8, K = 16, r = 4, t = 2, and Pf a = 10 -4 . 










































